Ultra high energy neutrinos from gamma ray bursts 
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Protons accelerated to high energies in the relativistic shocks that generate gamma ray bursts 
photoproduce pions, and then neutrinos in situ. I show that ultra high energy neutrinos (> 10 19 eV) 
are produced during the burst and the afterglow. A larger flux, also from bursts, is generated via 
photoproduction off CMBR photons in flight but is not correlated with currently observable bursts, 
appearing as a bright background. Adiabatic/synchrotron losses from protons/pions/muons are 
negligible. Temporal and directional coincidences with bursts detected by satellites can separate 
correlated neutrinos from the background. 



PACS numbers: 96.40 Tv, 
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The recent discovery of gamma ray bursts' (GRBs) af- 
terglows [Q, accurately predicted by theoretical models 
M and disappearance of flares in the radio flux || have 
bolstered our confidence in the correctness of the fireball 
model (J]. According to the model, bursts are generated 
when two or more hyperrelativistic shells, issued by an 
as yet unspecified source, collide with each other. A rela- 
tivistic shock forms, where non-thermal electrons are ac- 
celerated and then dissipate their internal energy through 
synchrotron (and possibly Inverse Compton) radiation. 
After the internal collision, the resulting shell will collide 
with the interstellar medium (ISM), thereby forming a 
second, relativistic shock, which will continue to expand 
into the ISM even after the burst proper, thusly generat- 
ing the afterglow. 

The relativistic environment surrounding the above- 
mentioned shock is ideal for the acceleration of protons to 
high energies @. The highest energy that can be attained 

is a 
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Here the explosion energy is E = -E52IO 52 erg, the expan- 
sion Lorenz factor 77 = 772 10 2 , the beaming angle 6, and 
the ISM number density n = n\ cm~ 3 . Currently popu- 
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lar values inferred from afterglows are 
§ implying e max « 6 x 10 20 eV. 

When energetic protons interact with synchrotron pho- 
tons emitted by electrons, they can produce pions; the 
decay of charged pions then produces electron and muon 
neutrinos. In this Letter, I will consider only ultra high 
energy neutrinos (UHENs, > 10 19 eV) and neglect lower 
energy ones ||,§. 

Expected fluxes of Ultra High Energy Neutrinos. Let 
us consider a burst of duration T seconds; according 
to the fireball theory for external shocks this oc- 
curs at a distance r e = rfcT from the unspecified burst 
source, and, in the shell frame, the shell thickness is 
5r = rjcT . The total energy density in the shell frame 
is then U, 



rcf. ll I find the inverse of the timescale for photopion 
losses, t~ , multiplying by the time the proton spends in 
the shell, in the shell frame (= rfT), I find that, for a pro- 
ton of energy e p as seen by an outside observer, immersed 
in a radiation field with turnover frequency e 7 rj 1 MeV, 
beyond which the spectrum significantly steepens, the 
total probability for photopion production is 



/(°) = 0.03% 
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for proton 



10 15 eVr)l(\ 



energies exceeding | «i = 2 x 
2 y± MeV/e^). I have used here a typical lu- 
minosity for long-lasting bursts, such as those with the 
ISM are thought to be, and a typical long duration. 

Experiments such as AIRWATCH |l^,|ll| have appre- 
ciable detection efficiencies for neutrinos exceeding the 
threshold energy e u j ss 10 19 eV. Since neutrinos emit- 
ted through photopion processes typically carry away a 
fraction q w 0.05 of the proton energy (losses will be 
discussed later), I have to compute the energy release in 
protons with energies exceeding e; = e Uy i jq k 2 x 10 20 eV. 
The spectrum in high energy protons accelerated at rel- 
ativistic shocks is roughly oc e~ 2 , and defining the total 
energy released in ultra high energy cosmic rays (UHE- 
CRs, e > ei = 10 19 eV) as E v , I have that the whole 
energy in UHECRs which can emit detectable UHENs 
(i.e., e > e t = e v .ijq) is E v In 

r max ^ in ax /ei. Only 

a fraction 2qf^ of this ends up in UHENs. Thus the 
total energy emitted in UHENs is 

In e r , 



E v = 2qf^E, 



(3) 
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The total flux of UHENs can then be obtained by inte- 
grating the flux over all distances: 

tigrbEu cK 
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7rr 2 c?7 2 ; inserting this into Eq. 3 of 
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where e v = e u j \ne max /ei is the average neutrino energy 
from this process, and the delicate factor K, to be dis- 
cussed later on, takes account of such unknowns as the 
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GRBs' redshift and luminosity distributions, and the de- 
tails of the cosmological model. 

The dependence of these neutrino rates upon physi- 
cal factors of individual bursts, such as T),Lj, and e 7 is 
all contained within (Eq. ||), and will be omitted 
from now on for sake of conciseness. The key factor in 
the above equation is E = tigrbEu, the injection rate 
per unit volume of non-thermal proton energy, because 
the others either are known or enter logarithmically. It 
is known already that, under the hypothesis that GRBs 
emit about as much energy in the form of 7-band photons 
and UHECRs, the flux of UHECRs at Earth is repro- 
duced to within a factor of w 3 0,0- I will show later 
that UHECRs are accelerated within afterglows, which 
dominate the energy balance by about a factor of 10. 
Then, if the same rough equipartition between radiation 
and UHECRs holds during the afterglow, the total energy 
release required to explain the UHECR seen at Earth is 
correctly accounted for. 

That the equipartition argument yields a correct an- 
swer can be checked by considering that the observed 
burst rate (~ 30 yr~ x Gpc~ 3 ) times the observed energy 
release including afterglow (« 10 52 erg) yields an energy 
release rate, 3 x 10 44 erg yr~ x Mpc~ 3 , very close to that 
deduced |l3| without explicit reference to the nature of 
the sources of UHECRs: E = 4.5 x 10 44 erg yr~ x Mpc~ 3 
for the restricted range of proton energies 10 19 eV < e < 
10 21 eV. 

Thus, under the equipartition assumption I can use the 
energy release necessary to explain Earth observations as 
the energy released in UHECRs by GRBs; taking e\ = 
10 19 eV, and defining H a = /i50 km s" 1 Mpc" 1 , I obtain 
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The flux determined above is not the whole flux of 
UHENs from GRBs detectable at Earth. The reason is 
that all UHECRs eventually will emit UHENs by pho- 
toproduction with photons of the CMBR, the so-called 
Greisen-Zatsepin-Kuz'min effect This neutrino pro- 
duction will occur in flight, rather than in situ, with a 
typical mean free path of order lOMpc. As they cross 
this distance, UHECRs are slowed down in their progress 
toward Earth by the turbulent intergalactic magnetic 
field. While estimates of this delay are very uncertain 
because of our ignorance of both strength and correla- 
tion length of the field, they still all agree in putting it 
above 10 2 — 10 3 yr, i.e. in washing away any correlation 
with GRBs observed within our lifetimes. The total flux 
of background UHENs h^ 9 \ uncorrelated with observ- 
able GRBs, is thus 



n^a) = ^,7.3x lO-^Kh-'yr- 1 cm~ 2 
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hypothesis on the distribution of redshifts and luminosi- 
ties of GRBs. A detailed computation Q for idealized 
redshift distributions of standard candles, has been car- 
ried out. Comparison of Table 1 of ref [jl5| with the above 
equation shows that their computed values of K vary by 
a factor of 3 either side of the value I obtained. 

Afterglows. I show now that acceleration of protons 
to the highest energies does continue unabated through 
most of the afterglow. After the burst, the relativis- 
tic shell keeps plowing through the interstellar medium, 
sweeping up more matter and decelerating. The shell 
Lorenz factor scales as r\ = %An x E^tJ s , where f<j is 
the post-burst time in days neglecting redshift. For adi- 
abatic expansion s = 3/8 while s — 3/7 for radiative 
expansion |17|. The maximum energy of non-thermal 
protons (Eq. 1) decreases very slowly with time, as i -1 / 8 
or t -1 / 7 for adiabatic or radiative expansion, respectively. 
In particular, for the best values E 52 = 1 and 9 = 1/3, 
production of UHENs ceases (i.e., e max < e v j/q) for 
7/ < 3.3, corresponding to w 6 d after the burst, nearly 
independent of whether expansion is adiabatic or radia- 
tive. 

I also have to check that the probability of photopion 
production through the afterglow does not change by 
much from the value computed (Eq. |^) for the burst 
proper. This requires some discussion. From observa- 
tions fisfl we know that the instantaneous luminosity 
scales as t~ a , with a « 1.1. Also, we know from fire- 
ball theory that T cx r/rj 2 , and that r oc r\~ v ', where 
v = 2/3 for adiabatic or v = 1/3 for radiative expan- 



sion. So the factor L 7 ?] T 
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However, it 



The computation of the factor K requires an explicit 



is more difficult to establish the variation of the spectral 
break e 7 with time, which is not currently observed. It 
seems however that, given the general softening of radi- 
ation within the afterglow, it is unlikely to remain con- 
stant; a more likely hypothesis is that it decreases slowly 
with time. Phenomenologically, one may take e 7 oc rf. 
The limits within which q is expected to vary are easy 
to ascertain. On the one hand, q — would imply that 
the cut-off does not evolve, despite the shell slow-down. 
This is both unphysical, and contrary to some weak evi- 
dence that it may decrease within the burst proper. On 
the other hand, the synchrotron turn-on frequency (i.e., 
that beyond which synchrotron emits most of the en- 
ergy) scales as oc j 4 ; in the afterglow model, all emission 
is due to synchrotron processes. However, the very long 
lasting optical emission from GRB 970228 seems to im- 
ply a very extended synchrotron spectrum, so that q = 4 
may be considered an upper limit. Thus < q < 4. I 
then obtain e 7 oc r) q oc t~ qs . From Eq. || I then find 
/i- "* oc t z , with z = s(q + 2 — v) — a. Only taking a small 
value, 5 = 1, and then only for adiabatic expansion, do 
I find z < 0. Thus we see that overall, the probabil- 
ity /i- - 1 is unlikely to decrease: if anything, is likely 
to increase through the afterglow, so that our estimates 
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are, most likely, lower limits. Thus, by taking in the 
previous section fj^ « constant, I did not overestimate 
the neutrino fluxes. An interesting consequence of this 
is that the luminosity in UHENs scales approximately as 
L v = £7 cc t _1 , which means that equal logarithmic 
post-burst-time intervals are equally likely to contain an 
observable neutrino. 

Losses. Proton losses (synchrotron and photo- 
hadronic) were shown to be negligible in ref. ||: the 
proton energy is limited by the size of the shell. I have to 
consider however adiabatic and synchrotron losses by pi- 
ons and muons, which could considerably limit the high- 
est energies achieved by neutrinos. 

Adiabatic losses are significant whenever the particle 
lifetime 7*7-* in the shell frame exceeds the characteristic 
timescale on which the magnetic field decreases because 
of the shell expansion; here * indicates either pion or 
muon, and tv = 2.6 x 1CP 8 s and t m = 2.2 x f CP 6 s are 
their respective lifetimes in their rest frames. The lim- 
iting Lorenz factors are found when the two timescales 
match, i.e., when 7*7"* = 2B/B. Following ref. JljJ 1 
shall take B cx i? -2 , where R is the transverse dimen- 
sion of the causally connected region, which, following 
refs. ]f9| , p0| , is given by r/rj, even through the afterglow, 
and obviously R k, c. Then I obtain the limiting Lorenz 
factor in the observer frame 7/ = r/cr*, independent of 
whether the afterglow is adiabatic or radiative. Scaling 
r = xri by its lowest value, that at the moment of the 
burst proper, n = 2r] 2 cT = 6 x 10 15 cmr£(T/ 10 s), 1 
find 7^ = 10 13 ir, and 7^ = 10 n :r for x > 1, both ex- 
ceeding the proton's Lorenz factor in Eq. 1. For protons 
with Lorenz factor j p in the shell frame, the synchrotron 
cooling time is t s = 1 yr (10 n /7p)(l G/B) 2 . For syn- 
chrotron losses to be negligible, the Lorenz factor of pi- 
ons/muons must not exceed the limiting 7* given by |lE(] 
'ypt s (m+/m p ) 3 = 7 2 r*, where m+/m p ss 0.1 for both pi- 

ons and muons. From ref. p0[ |, B » 1 G r\^ 2 for the 
external shock scenario and the afterglow. Transforming 
back to the observer frame 1 find 7„ = 3x 10 13 ??2 /2 , and 
7/i = 3 x lO 12 ??^ 2 . Both exceed the Lorenz factors of 
the proton, Eq. 1. Thus adiabatic/synchrotron losses of 
pions/muons do not affect the arguments of this paper. 

Detectability. Currently planned experiments such as 
A1RWATCH will monitor from satellites fluorescent 
light profiles of cosmic ray cascades over areas of order 
A = As x 10 6 km 2 , with A& k 1. The interaction proba- 
bility for UHENs is proportional to the monitored column 
density (10 3 g cm~ 2 ); it also depends over the extrapola- 
tion of the cross-section to currently unobserved energies, 
but typical values are a w 3 x 1(T 32 cto 2 (e^/10 19 eU) 1 / 2 
jll],|l]]. Once the neutrino has interacted, a detection 
efficiency close to 1 for UHENs is reported by feasibility 
studies, at energies e„ w 10 19 eV, yielding interaction 
probabilities of P v m 3 x 10~ 5 . This translates into an 
expected number of detectable UHENs of 



AO) 

b y 0.03 
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At the same time, we expect a background flux from Eq. 
I given by 



N( b9) = 200 K A 6 yr^h- 1 



(8) 



It is safe to state that Eqs. ^ and ^ have large errors, 
due to our ignorance both of the neutrino-nucleon cross- 
section at these large, and untested neutrino energies, 
and to the sources' redshift distribution (the parameter 
K). 

The requirement that the expected number of neutri- 
nos correlated with bursts be large enough to ensure de- 
tection within a year of operation can be turned, using 
Eq. ||, into a requirement on the area covered by the 
experiment: 
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Detection of correlated neutrinos seems possible provided 
bursts due external shocks are well-represented by the 
average values employed above. 

The flux of Eq. |^ of an event per day, completely 
uncorrelated with currently observable bursts, obliges us 
to face the issue whether we can distinguish from casual 
associations a much smaller (/i° ~ 0.03) flux which is in- 
deed correlated (to within the afterglow duration, as 6 d) 
with simultaneously observed bursts. The answer would 
be an easy yes if UHENs arrived simultaneously with the 
burst proper, because we could then use very tight direc- 
tional and temporal coincidences to distinguish the signal 
from background noise. But, since I argued above that 
most neutrinos are produced during the afterglow which 
is observed to last for a few days after the burst, it has 
to be ascertained whether this can still be done. The 
answer is a qualified yes. 

Suppose I can measure the directions of arrival of neu- 
trinos and GRBs with a combined directional error of [3. 
Calling Ngrb the rate of detection of GRBs in the 7-ray, 
the probability of casual association P c is 



2 

N GRB db h — = 4 x 10~ 



6 d V 1 



(10) 



where I used N GRB = 300 y^ 1 , typical of BATSE 
p2| . The rate of appearance of casual associations is 
N^ 3 ^ P c w 0.08 yr~ x , reassuringly smaller than the 
rate of physical associations, Eq. 0. This condition, 
PcN^ 9 ^ N„, can also be written as 



St [13 



/<°> » Ngrb&^jm 4xl0- 4 f- d \f o 
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Comparison with Eq. g shows that the experiment can 
be done, provided angular errors of order 
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3<£-7°( ^ h — — \ " (12) 

' V 10 50 erg rfo 1 MeV 10 s 6d J K ' 

can be achieved. 

Lastly, since the rate of Eg. p| is comparable to that of 
GRBs detected by BATSE [ p2[ , measurement of dipole 
and quadrupole moments of the neutrino distribution 
may just be doable. The spectrum of UHENs (both back- 
ground and correlated ones) will follow accurately that 
of UHECRs in GRBs, since the probability of photopion 
losses (Eq. ||) is independent of proton energy. It should 
thus be possible to see the cutoff in the UHECR spec- 
trum, Eq. 1, as mirrored in neutrinos. 

Discussion. The acceleration of UHECRs in GRBs 
is so effective, that it has been proposed |^,|l2| that the 
whole flux of UHECRs at Earth comes from these events. 
However, since UHECRs can take ~ 10 3 yr more than 
photons to reach us from the closest GRBs, it will be 
impossible, within our finite lifespans, to establish a di- 
rect association between GRBs and UHECRs. A sure 
hint should be that no AGNs, or peculiar object, ought 
to be seen close to the direction of arrivals of UHE- 
CRs, but this expectation is not unique to this model, 
and is common for instance to strings. On the other 
hand, a UHEN of rj 10 19 eV would accumulate with re- 
spect to photons emitted simultaneously a delay of only 
rs 1CU 19 s(m u /l0 eV) in coming from even a distance 
of c/H , the radius of the Universe, with m v the neu- 
trino mass. Thus it would be essentially simultaneous to 
photons (including afterglow's photons). Furthermore, 
the UHENs can only be produced by the highest energy 
protons, those, in other words, well beyond the Greisen- 
Zat'sepin-Kuzmin limit. Thus the UHENs produced in 
situ represent the surest smoking gun that UHECRs are 
accelerated in GRBs. Different, electromagnetic signa- 
tures of the association of UHECRs and GRBs have been 
discussed in refs. |23| , p^| . 

In short, what detection of UHENs will allow us to do 
is to circumvent the shortsightedness imposed upon us by 
the Greisen-Zatsepin-Kuz'min limit, and to investigate 
the generation of the highest energy cosmic rays through- 
out the whole Universe. The only alternative sources of 
UHENs proposed so far are cosmic strings [^5| and AGNs 
p6| which are also the only alternative sources proposed 
so far for UHECRs. I have discussed here that a fraction 
of all UHENs (Eq. ^) should show an association with 
simultaneously observed GRBs, if they indeed originate 
in GRBs. Thus a potentially clear-cut way to distinguish 
between the three competing theories is available and it 
might, perhaps, already be accessible to AIRWATCH- 
class experiments. 

Thanks are due to L. Scarsi and J. Rachen for helpful 
comments. 



[1] In the X-ray, Costa, E. et ai, 1997, Nature, 387, 783; in 
the optical, van Paradijs, J., et al, 1997, Nature, 386, 
686; in the radio, Frail, D.A., et ai, 1997, Nature, 389, 
261. 

[2] In the radio, B. Paczyhski, and J.E. Rhoads, Ap.J.L., 

418, L5 (1993); in the optical, P. Meszaros, and M.J. 

Rees, Ap.J., 476, 232 (1997); in the X-ray, M. Vietri, 

Ap.J.L., 478, L9 (1997). 
[3] J. Goodman, New Astronomy, 2, 449 (1997). 
[4] For reviews, see T. Piran, in Some unsolved problems in 

Astrophysics, Princeton Univ. Press, Princeton 1995, eds. 

J. Bahcall and J.P.Ostriker, p. 435; M.J. Rees, astro-ph 

n. 9701162 (1997). 
[5] J.J. Quenby, Lieu, R., Nature, 342, 654 (1989). 
[6] M. Vietri, Ap.J. 453, 883 (1995); the dependence of Eq. 

1 upon n\ differs from that in this reference because a 

small error has been corrected. 
[7] E. Waxman, S. Kulkarni, D.A. Frail, submitted to 

Ap.J.L., astro-ph. n. 9709199 (1997). 
[8] B. Paczyhski and G. Xu, Ap. J., 427, 708 (1993). 
[9] E. Waxman and J. Bahcall, Phys. Rev. Lett., 78, 2292 

(1997). 

[10] J. Linsley, in Proc. of the XXVth ICRC, Durban RSA, 5, 
381 (1997). 

[11] Y. Takahashi, et ai, SPIE 2806-13 (1996); J. Linsley et 

ai, in Proc. of the XXVth ICRC, Durban RSA, 5, 385 

(1997); P. Attina et at, ibidem, 5, 389 (1997); J. Forbes 

et al, ibidem, 5, 273 (1997). 
[12] E. Waxman, Phys. Rev. Lett., 75, 386 (1995). 
[13] E. Waxman, Ap. J. L. 452, LI (1995). 
[14] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G.T. Zat- 

sepin and V.A. Kuz'min, JETP Lett. 4, 78 (1966); R.J. 

Protheroe and P.A. Johnson , Astroparticle Phys. 4, 253 

(1995) and erratum 5, 215 (1996). 
[15] S. Yoshida and M. Teshima, Prog. Theor. Phys., 89, 833 

(1993). 

[16] E. Maoz, 1994, ApJ, 428, 454 (1994). 

[17] M. Vietri, Ap.J.L., 488, L105 (1997). 

[18] A. Fruchter et al, IAU Circ. n. 6747 (1997). 

[19] J. P. Rachen, and P. Meszaros, submitted to Phys. Rev. 

D, astro-ph. n. 9802280 (1998). 
[20] P. Meszaros, P. Laguna, and M.J. Rees, ApJ, 405, 278 

(1993). 

[21] C. Quigg, M.H. Reno, and P.T. Walker, Phys. Rev. Lett., 
57, 1204 (1986). 

[22] G.J. Fishman and C.A. Meegan, ARAA, 33, 415 (1995). 

[23] M. Vietri, Phys. Rev. Lett., 78, 1323 (1997). 

[24] M. Boettcher and CD. Dermer, sumitted to ApJ, astro- 
ph. n. 9801027 (1998). 

[25] G. Sigl, S. Lee, D. Schramm and P. Bhattacharjee, Sci- 
ence, 270, 1977 (1995). 

[26] J. P. Rachen and PL. Biermann, Astron. Ap., 272, 161 
(1993). 



4 



